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INTRODUCTION 


One  of  the  questions  that  arises  in  the  area  of  straightening  components 
deals  with  the  magnitude  of  deflections  and  the  corresponding  strains  which 
could  adversely  affect  material  properties  and  limit  subsequent  usage  of  the 
component.  An  equally  pertinent  question  centers  on  the  need  for  hot 
straightening  to  eliminate  excessive  permanent  deflections  or  bows  in  compo¬ 
nents  due  to  processing  operations.  With  these  questions  forming  a  backdrop, 
the  writer  performed  an  analysis  comparing  three  and  four  point  bending  from 
the  standpoint  of  seeing  whether  larger  permanent  deflections  could  be  removed 
using  the  four  point  method  and  thus  possibly  obviating  the  need  for  hot 
straightening.  It  was  hoped  that  this  could  be  done  without  causing  addi¬ 
tional  material  degradation  as  a  result  of  straining  the  material  into  the 
plastic  region. 

The  intent  of  this  report  is  to  present  the  results  of  a  theoretical 
elastic-plastic  analysis  as  a  means  of  comparing  the  two  methods. 

The  theory  used  as  a  basis  for  the  calculations  in  this  report  is 
essentially  the  same  as  developed  by  Seely  and  Smith. 1  We  first  consider  the 
case  for  a  concentrated  load  at  midspan  as  shown  in  Figure  1(a).  Using  the 
dummy  load  method,  which  is  practically  the  same  as  the  principle  of  virtual 
work,  the  deflection  can  be  calculated  from  the  integral 

6  =  J  mda  (1) 


If.  B.  Seely  and  J.  0.  Smith,  Advanced  Mechanics  of  Materials,  2nd  Ed.,  Wiley 
and  Sons,  1952. 
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where  m  is  the  moment  due  to  the  dummy  load  and  da  the  relative  rotation  of 


one  plane  of  a  cross  section  with  respect  to  another  a  distance  dx  apart.  For 
the  linear  region  da  is  equal  to  (M/EI)  dx,  where  M  is  the  moment  due  to  the 
real  load,  E  is  the  modulus  of  elasticity,  and  I  is  the  moment  of  Inertia.  We 
first  integrate  through  the  elastic  region  from  0  to  X,  where  X  is  the  dis¬ 
tance  from  the  left  support  to  the  end  of  the  elastic  region.  We  then  have 


(2) 


The  limit  X  depends  on  the  loading  conditions  and  can  be  obtained  by  ratio 
from  the  moment  diagram.  Referring  to  Figure  1(b) 
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FTgura  Load  and  Moment  Diagrams  for  Three  Point  Bending. 
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As  the  moment  increases  from  the  yield  moment  My  to  the  fully  plastic  moment 
Mfp,  X  decreases,  indicating  that  the  length  of  the  elastic-plastic  zone 
increases . 

We  next  integrate  through  the  elastic-plastic  zone  from  X  to  i/2,  For 

the  elastic-plastic  region  da  now  takes  the  form  of 

1  /  2  . 

Ay  n  dx 


da  = 


EI  /K  -  Mp 


(5) 


Mv 


One  can  refer  to  the  Appendix  for  the  development  of  this  expression.  Hence, 


62  =  / 


111  %  n^^mdx 

x  EI  /K  _  Mp 


(6) 


M„ 


The  total  deflection  can  then  be  obtained  by  adding  deflections  6^  and  62  and 
multiplying  by  two  since  we  have  symmetrical  loading. 

For  a  simply  supported  beam  having  a  concentrated  load,  M  -  Px/2  and  m  = 
x/2  for  the  region  0  <  x  <  i/2 •  Hence: 


X  P  x2 

61  =  /  -  —  dx 

0  4  EI 


%  ,  /9  Ml  x 

and  6o  =  —  n1'2/  - 

EI  X  2 


/K  -  Px 


dx 


(7) 


(8) 


2MV 


Inserting  the  factor  two  previously  mentioned,  we  obtain: 

x 

vV  _  Px 


P  .X  „  My  ,  ,n  rHl 


6  =  6i  +  6?  = - I  x2dx  +  —  n-*-/2  / 

2EI  0  EI  X 


dx 


(9) 


2M,; 


3 


Integrating  and  inserting  limits,  we  get  the  equivalent  expression  given  by 
Seely  and  Smith! 


6  “  TTZ  t1  +  2n(2K+l)  -  2nl/2(2K  +  ~->(K  -  1 


PA 


1/2 


48EI  Mniax 


4Hv 


4Mv 


(10) 


It  is  readily  apparent  that  this  expression  degenerates  to  the  elastic 
equation  when  M^x  =  My.  That  is, 


6 


Pi.3 

48EI 


(ID 


A  computer  program  was  written  for  equation  (10)  and  given  the  ancronym  DEFL 
The  program  has  a  Do  Loop  on  P.  The  results  shown  graphically  were  obtained 

from  this  program. 

We  next  consider  the  case  for  two  concentrated  loads  symmetrically 
located  relative  to  the  midspan  as  shown  in  Figure  2(a). 


Figure  2.  Load  and  Moment  Diagrams  for  Four  Point  Bending. 


Ip .  B.  Seely  and  J.  0.  Smith,  Advanced  Mechanics  of  Materials,  2nd  Ed.,  Wiley 
and  Sons,  1952. 
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In  this  case,  we  have  to  integrate  through  three  regions,  0  to  A,  A  to  A,  and 
A  to  i/ 2. 


X  Mmdx  1  A  x 

61  =  J  - =  —  /  Px(-)dx 

0  El  El  J0  2 


(12) 


Mv 


where  the  limit  X  by  ratio  is  equal  to  -  A  as  can  be  noted  from  Figure 


M, 


2(b). 


■max 


A  My  n1//2mdx  My  .  A  x/2 
62  =  J,  ~  - -  —  n1'2  /  - dx 


X  EI  v^k  -  Mp  EI  *  /K  _  Px 


(13) 


M-u 


and 


Mt, 


i/2  n^^mdx 


63  =  / 


EI  /K  -  ”®ax  EI  /—  H,ax  A 


,i/2  x 
/  (-)dx 


(14) 


M,; 


M-u 


Evaluating  these  integrals,  and  inserting  the  factor  of  two  because  of 

symmetry  and  integrating  over  half  the  length,  the  maximum  deflection  at  i/2 
becomes 


6  —  6^  +  $2  +  (S3 


(15) 


1/2 


1/2 


6  = 


PA3  My3  2  My3  „  tn 

55  W  1  ir  ir  +  JT>«  -  “>  -  (2K+l)(K-l)l/2] 


My  n1/2  rl2 


+  —  “ZZZZZZ.  [ 
EI  /K  -  pa  8 

M„ 


(16) 


A  computer  program  was  written  for  equation  (16)  and  given  the  ancronym 
SSDEFL2P.  The  program  has  a  Do  Loop  on  P.  The  results  shown  graphically  were 
obtained  from  this  program. 
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If  we  specialize  this  equation  for  A  =  A/2,  the  second  and  third  terms 
drop  out  and  we  get 

6  =  1-1  (17) 

24EI 

This  differs  from  equation  (11)  by  a  factor  of  two  since  we  have  two  loads 
acting  at  A/2. 

RESULTS  AND  DISCUSSIONS 

As  a  means  of  comparing  three  point  bending  with  four  point  bending,  the 
following  cases  were  calculated  using  the  previously  mentioned  computer 

programs  with  the  following  input  data: 

Ro  =  3.0,  3.5,  4.0,  4.5  inches 
=  2.0  inches 
E  =*  30  x  106  psi 
oy  =  160  Ksi 
A  *  60  inches 
A  =  15,  20,  25  inches 

The  load  P  was  located  at  A/2  for  the  case  of  three  point  bending.  This  gave 
a  total  of  16  load-deflection  curves.  An  example  for  the  case  of  R0  =  3.0 
inches  with  A  *  15,  20,  and  25  inches  is  shown  in  Figure  3.  Assuming  linear 
elastic  unloading,  the  permanent  deflections  were  obtained  by  projecting  down¬ 
ward  from  the  maximum  deflection  at  maximum  load,  to  the  deflection  axis  with 
a  line  parallel  to  the  initial  elastic  loading  line.  It  should  be  pointed  out 
that  load  P  which  would  be  applied  to  the  loading  fixture  to  accomplish  the 
four  point  loading  is  actually  equal  to  2P  and  is  consequently  twice  the  value 
of  P  applied  in  the  three  point  case.  The  same  is  also  true  for  all 
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subsequent  diagrams  and  figures  comparing  three  point  and  four  point  loading. 
From  these  16  curves  we  constructed  two  others,  Figures  4  and  5*  Figure  4  is 
a  plot  of  moment  of  inertia  (I)  vs.  deflection  for  the  three  point  loading  and 
three  different  values  of  A  for  the  four  point  loading.  The  values  of  I  for 
four  different  outside  radii  are  marked.  The  solid  lines  show  the  maximum 
deflection  under  load  and  the  dashed  lines  are  for  the  permament  deflections. 
The  plasticity  condition  is  such  that  the  elastic-plastic  interface  has  moved 
a  distance  equal  to  75  percent  of  the  outside  radius  from  the  outside  fiber 
toward  the  neutral  axis.  Equation  35  of  reference  2  gives  a  relationship 
between  the  outside  radius  Rq,  the  distance  from  the  neutral  axis  to  the 
elastic-plastic  interface  p,  and  the  yield  strain,  ey.  That  is, 

Ro 

emax  ~  "  Gy 

R0~P 

For  the  75  percent  plasticity  condition  used  for  comparison  purposes  for  the 
two  methods,  Rq,  p,  and  Gy  would  be  the  same,  therefore  the  maximum  strains 
would  be  the  same  for  both  bending  methods.  The  same  reasoning  would  hold  for 
the  nearly  100  percent  plasticity  condition  also  considered. 

Figure  5  is  the  same  plot  as  Figure  4  except  the  plasticity  condition  is 
such  that  the  elastic-plastic  interface  has  moved  a  distance  equal  to  nearly 
100  percent  of  the  outside  radius.  Table  I  gives  a  summary  of  the  various 
geometrical  parameters,  loads,  moments,  and  deflections. 


2r.  V.  Milligan,  " Moment-Strain  Relationships  in  Elastic-Plastic  Bending  of 
Beams,"  to  be  published. 
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Figure  6  was  obtained  by  taking  ratios  of  maximum  deflections  under  load 
for  the  four  point  case  to  the  three  point  case  for  a  constant  R.  For  the  100 
percent  plasticity  condition,  the  maximum  ratio  is  7.85  for  A  =*  15  and  Rq  ■ 

3.0  while  the  minimum  ratio  is  3.25  for  A  =  25  and  R 0  38  4.5  .  For  the  75 
percent  plasticity  condition,  the  maximum  ratio  is  2.90  for  A  *  15  and  R0  50 
3.0  while  the  minimum  ratio  is  1.61  for  A  3  25  and  Rq  *  4.5.  It  is  thus 
evident  for  both  plasticity  conditions  that  the  maximum  ratios  occur  for  the 
beam  having  a  smaller  cross  section  and  as  the  loads  get  nearer  to  the 
supports • 

Figure  7  is  the  same  as  Figure  6,  except  the  deflection  ratios  are  perma¬ 
nent  deflections.  Again  the  maximum  ratio  for  the  100  percent  plasticity  con¬ 
dition  occurs  for  A  -  15  and  Rq  *  3.0  while  the  minimum  ratio  occurs  for  A  * 

25  and  Rq  =*  4.0.  For  the  75  percent  plasticity  condition,  the  maximum  ratio 
occurs  for  A  *  15  at  a  value  of  R0  between  3.5  and  4.0. 

In  summary,  the  ratio  of  maximum  deflections  for  the  100  percent 
plasticity  condition  run  from  a  low  of  about  3.25  to  a  high  of  about  7.75. 

For  the  75  percent  plasticity  condition  they  run  from  a  low  of  about  1.75  to  a 
high  of  about  3.75.  The  ratios  are  even  greater  for  the  permanent 
deflections.  Here  they  run  from  a  low  of  6.25  to  a  high  of  19  for  the  100 
percent  condition  and  from  about  2  to  8.75  for  the  75  percent  condition.  It 
is  thus  very  evident,  based  on  the  theory  presented  here,  that  four  point 
bending  can  be  a  means  of  significantly  increasing  the  deflections  compared  to 
three  point  bending,  while  at  the  same  time  keeping  the  strains  and  material 
degradation  at  the  same  level. 


In  an  effort  to  propose  a  reason  for  the  apparent  advantage  of  four  point 
bending,  the  author  submits  the  concept  of  distributed  plastic  flow.  Figure  8 
is  a  plot  of  maximum  deflection  under  load  vs.  length  of  plastic  zone.  The 
length  of  plastic  zone  can  be  determined  by  simple  ratio  as  follows.  For  the 
four  point  case: 

X  A  My  My 

—  -  - -  hence  x  - - -  A  =  - - 

%  ^max  pmaxA  pmax 

where  s  is  the  distance  from  the  left  support  to  the  beginning  of  the  plastic 
zone.  For  the  three  point  case: 

X  i/2  My  My 

—  = -  hence  x  - - (i/2)  - - 

My  Mjnax  Mnjax  Pfnax/ 2 

Figure  8  shows  sketches  depicting  the  size  of  the  elastic-plastic  zones 
both  for  the  three  point  case  and  three  different  values  of  A  for  the  four 
point  case.  Obviously,  as  the  loads  move  outward  toward  the  supports,  the 
elastic-plastic  region  becomes  larger.  The  plot  of  vs.  length  of  plastic 

zone  is  nearly  linear.  This  tends  to  support  the  concept  that  as  one  distrib¬ 
utes  the  plastic  flow  one  can  increase  the  deflection  without  appreciably 
increasing  the  maximum  fiber  strain  with  its  probable  consequence  of  degrada¬ 
tion. 

CONCLUSIONS 

From  the  analysis  presented  it  appears  that  larger  deflections  under  load 
can  be  obtained  and  consequently  larger  permanent  deformations  can  be  removed 
by  using  the  four  point  bending  method  compared  with  three  point  bending  while 
keeping  the  maximum  fiber  strains  the  same  in  both  cases. 
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The  concept 
the  advantage  of 


of  distributed  plastic  flow  appears  to  be  a  valid  reason  for 
the  four  point  method  relative  to  the  three  point  method. 
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APPENDIX  A 


In  a  previous  report^  the  bending  moment  corresponding  to  an  elastic- 
plastic  interface  depth  p  for  a  rectangular  beam  was  developed  and  shown  to 
be 


Oyb 

Mp  =  “7“  [h2  +  2ph  -  2p2]  (Al) 


Figure  Al .  Strain  vs.  Depth  of  Elastic-Plastic  Interface. 
From  Figure  Al  we  get  the  following  relationship 

ep  1 


(  A2) 


Solving  for  p ,  we  obtain 


1  -  2p 

h 


h  ey 

P  =  -  (1 - ) 

2  ep 


(A3) 


Now  substituting  this  into  equation  (Al) ,  we  get  the  bending  moment  correspon¬ 
ding  to  the  elastic-plastic  interface  being  at  a  depth  equal  to  p. 


Oybh2  3  1  eyZ 

mp  -  —  b  - ;  ri\ 


2  £r 


(A4) 


^R.  V.  Milligan,  "Moment-Strain  Relationships  in  Elastic-Plastic  Bending  of 
Beams,”  to  be  published. 
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Dividing  this  expression  by  the  yield  moment  My  =  bh2/6  we  get  the  moment 


ratio 


Mn 


M 


3  1  ey 

2  2 


( A5) 


Multiplying  and  dividing  the  second  terra  of  equation  (A5)  by  (dx/-)  we  get 


j'  Cydx  ^ 

Mp  3  i  h/2 

My  2  2  f £pdx^ 2 

MVT 

Now  using  the  results  of  equations  (A4)  and  (A 5)  we  have 


(A6) 


M 


p  3  1  (dOty)' 


(A7) 


My  2  2  (da)2 

With  the  view  of  attempting  to  make  equation  (A7)  more  general  so  that  it 
can  be  used  for  beam  cross  sections  other  than  rectangular,  Seely  and  Smith 
substitute  K  for  3/2  and  n  =  K  -  1  for  1/2.  K  now  represents  a  geometrical 
shape  factor  equal  to  Mfp/My,  i.e.,  the  ratio  of  the  fully  plastic  moment  to 
the  yield  moment.  Table  6.1  of  Faupelfs  book^  gives  values  of  K  for  different 
beam  cross  sections.  Substituting  these  parameters,  we  have 


Mf)  ^day^9 

--  -  K  -  n(-~)2 


M 


da 


(A8) 


3j.  H.  Faupel,  Engineering  Design,  Wiley  and  Sons,  1964,  p.  395. 
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Solving  for  da 


rd “y12 

K - -  n( - 2 

My  da 

day2 

da2  =  n  - 

Mv 


da  = 


n1/2  0»y 


( A9) 


Using  Hooke’s  law  and  the  flexure  formula  and  substituting  for  e  in  (B4)  we 
have 

My 

dav  =  —  dx 
y  El 


Substituting  this  into  equation  (A9)  we  obtain  the  desired  expression: 


da  = 


%  n1/2_dx 
EI 

Mv 


(A10) 
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APPENDIX  B 


At  incipient  yielding  on  the  outside  fiber  we  can  determine  an 
expression  for  da  as  follows.  By  referring  to  Figure  B1  which  indicates  that 
planes  remain  plane  and  further  assuming  that  the  tangent  of  the  angle  equals  ' 
the  angle  for  small  deformations,  we  have 


Figure  B1 .  Strain  vs.  Depth  of  Cross-Section. 

<5y 

tan(day)  =  dctv  - - 

7  7  h/2 

From  an  expression  for  engineering  strain,  we  obtain 

(dx+fiy)  _  dx  6y 

ey  =  =  — 

dx  dx 


hence 

substituting  into  equation  (Al),  we  get 

edx 


(Bl) 


(B2) 


(B5) 
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